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Abstract 



A method for computing low-temperature series for renormalized operators in the 
two-dimensional Ising model is proposed. These series are applied to the study of 
the properties of the truncated renormalized Hamiltonians when we start at very low 
^ ^ temperature and zero field. The truncated Hamiltonians for majority rule, Kadanoff 

transformation and decimation for 2x2 blocks depend on the how we approach the 
first-order phase-transition line. These Renormalization Group transformations are 
multi-valued and discontinuous at this first-order transition line when restricted to 
some finite-dimensional interaction space. 
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1 Introduction 



The behavior of the Renormalization Group (RG) in the vicinity of first-order phase 
transitions has been a very controversial matter for the last 20 years. In 1975 Nienhuis 
and Nauenberg [1] proposed that the RG transformations (RGTs) behave near first- 
order transition points in a similar fashion as near standard critical points. Each RG step 
is smooth (i.e. the renormalized couplings are analytic functions of the original ones, 
even at the transition points). Singular behavior is recovered as we infinitely iterate this 
transformation near the attracting manifold of a fixed point. Moreover, the fixed point 
governing a first-order phase transition has the following properties 1 : i) A domain of 
attraction which includes the first-order transition surface, ii) Zero correlation length at 
the fixed point (In most systems, first-order transition points possess a finite correlation 
length. See [2] for a counterexample), iii) A relevant operator whose critical exponent 
is given by the dimensionality of the system, y = d. As a matter of fact, there are as 
many exponents y = d as phases coexist at the transition line 2 [3]. In the Ising model 
it is believed that this fixed point is located at zero temperature [4]. 

This picture was criticized by some authors [5,6,7,8,9] who claimed that the RG flow 
is itself discontinuous at the transition line. That is, they claimed that the renormalized 
Hamiltonian has different limiting values depending on how the original Hamiltonian 
approaches the transition line. Most of these claims were based on Monte Carlo Renor- 
malization Group (MCRG) computations [5,6,7,8]. In ref. [9] non-rigorous analytical 
arguments were given to support the same conclusion. 

In opposition to this view, it was argued in ref. [10] that the observed discontinuities 
are due to the truncation of the Hamiltonian space inherent in the MCRG approach. 
That is, the exact RG map is continuous at the first-order phase transition (in agreement 
with the conventional scenario), but truncation could induce the observed discontinu- 
ities. In fact, for the two-dimensional (2D) Ising model and majority rule with 2x2 
blocks it was found that the discontinuity in the magnetic field was of the same order as 
the truncation error. Moreover, as the number of operators included in the computation 
was increased, the size of this discontinuity decreased. 

This puzzle was partially solved by van Enter-Fernandez-Sokal [11], who showed that 
for systems with bounded dynamical variables and interacting through a Hamiltonian 
belonging to the space B 1 (i.e. the space of real, absolutely summable and translation- 
invariant interactions) the RG map is always continuous and single- valued, whenever 
it exists at all (subject to some very mild locality conditions on the RGT). For finite 
systems the existence of the transformation (i.e. of the renormalized Hamiltonian) is 
trivial. In the thermodynamic limit, however, this is a very subtle problem. As a matter 
of fact, van Enter et al proved that the renormalized Hamiltonian does not exist in the 
2D Ising model when the temperature is low enough, for the Kadanoff transformation, 
decimation, block average and some particular cases of majority rule. The possibility 
of a non-Gibbsian renormalized measure had been explained earlier by Griffiths and 

1 Nienhuis and Nauenberg called this a discontinuity fixed point (DFP). However, we shall avoid this 
name as it causes confusion with a different scenario (discontinuous RG flow) described in the next 
paragraph. 

2 Here we take into account the (trivial) critical exponent y = d associated with the renormalization 
of the identity operator in the Hamiltonian. 
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Pearce [12] using more physical (but non-rigorous) arguments. On the other hand, 
what happens for the majority rule with 2x2 blocks (the case most considered in the 
literature) is still unclear. 

It therefore remains to be understood (a) why MCRG methods show a discontinuous 
RG flow at or near first-order phase transition points and (b) whether such apparent 
discontinuities are in any way connected with the Gibbsian or non-Gibbsian nature 
of the renormalized measure. The logic of the MCRG approach is the following: 1) 
Compute expectation values of suitable renormalized local operators using a Monte 
Carlo (MC) algorithm and a certain real-space RGT. 2) Assume that the renormalized 
measure is a Gibbs measure for some (a priori unknown) interaction 7i' n 3 belonging to a 
(pre-chosen) n-dimensional subspace V n C B 1 . 3) Compute the renormalized coupling 
constants (i.e. compute 7i' n ) using some statistical inversion method. In general, one 
has to minimize some strictly convex functional F n 4 , which leads to a highly non-linear 
set of equations involving the expectation values computed in step 1. MCRG methods 
thus study a truncated RG map in which the renormalized Hamiltonian 7i' n is forced 
to lie in the subspace V n . At this point we can ask two different questions: (A) For 
each fixed n < oo, what are the properties of 7i' n : existence, uniqueness, continuity,. . . ? 
(B) What is the behavior of 7i' n as n tends to infinity? Presumably, the answers to 
(A) and (B) are model-dependent: they depend generically on the physical system, on 
the RGT and on the choice of the subspaces V n . On the other hand, the connection 
between Gibbsianness and the behavior of the 7i' n as n — > oo is far from clear. In [11] 
it was proven that if the renormalized measure is Gibbsian, then there is a sequence of 
truncated interactions {7i' n } } which are almost minimizers of F n and which converge to 
the true renormalized Hamiltonian Ti 1 in B 1 norm. However, it is not guaranteed that 
the exact minimizers 7i' n also converge to 7i' . On the other hand, if the renormalized 
measure is non-Gibbsian, then it was proven in [11] that the the sequence {Ti-' n } has no 
limit at all in the space B 1 . 

The simplest system which undergoes a first-order phase transition is the 2D Ising 
model at T < T c . In particular, for T <C T c we can do a complete analytical study using 
low-temperature (low-T) expansions [13] 5 . Analytical methods would be preferred in 
this kind of investigation, as all sources of errors are under control 6,7 . 

In this paper we address question (A) for the 2D Ising model at low temperature, 
using series expansions. We develop a simple procedure to compute these series for the 
expectation values of the renormalized operators which enter in the equations needed to 
obtain 7i' n . For real-space RGTs the expectation value of an operator with respect the 

3 We shall hereafter denote renormalized quantities with a prime. 

4 In ref. [11], Section 5.1.2, the interaction %' n comes from minimizing the relative entropy density with 
respect to the renormalized measure. See also Section 3 below. 

5 It ought to be also possible to study the first-order transition undergone by the (/-state Potts model 
in the limit q 1 using l/q expansions. 

6 Series expansions do not suffer from two types of errors inherent to any MC simulation: statistical 
errors and finite-size corrections. These errors affect the estimators of the expectation values obtained 
in the MC simulation and they are propagated to the interaction %' n . By contrast, the predictions for 
such expectation values given by series expansions are obtained directly in the thermodynamic limit and 
no stochastic process is involved. 

7 Note that unlike many applications of series expansions, here we are really interested in the low-T 
behavior and not in the critical region T ps T c . Therefore, no extrapolation procedure is involved. 
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renormalized measure can be written as an expectation value in the original measure of a 
certain composite operator 0. This is equal to the original operator acted upon by 
a probability kernel (which is the mathematical object representing the RGT). Thus, if 
we know how to obtain the low-T expansions in the original (unrenormalized) measure, 
then we can compute any expectation value by doing the corresponding integral. 

These series can be useful in two other ways. They provide a check for MCRG com- 
putations at low temperature, since expectation values coming from the MC simulations 
can be compared with the low-T predictions. On the other hand, when performing a 
RGT the system is viewed at a larger spatial scale. For that reason we believe that the 
low-T series for the renormalized magnetization, susceptibility and specific heat could be 
used to extract the critical exponents (using standard series-extrapolation techniques). 
In fact, a better convergence could be expected for these "improved" series. It would be 
interesting to devise a computational procedure to generate these series to an arbitrary 
order. 

The goal of this paper is thus to study the properties of the finite-dimensional ap- 
proximants 7i' n for the 2D Ising model. Starting on the first-order transition line (i.e. 
at zero magnetic field) and at very low temperature, for each of the two pure phases 
we obtain (via some RGT) two different renormalized measures v'^\ For each 
of them we can find a unique truncated renormalized Hamiltonian 'H'^ ■ The natural 
question is: are these two truncated Hamiltonians equal or not? Or equivalently: do 
all the odd couplings in 1~C^ vanish, or not? An affirmative answer implies that the 
truncated RGT restricted to V n is continuous and single- valued at the transition line 8 . 

We have studied three different RGTs: decimation, Kadanoff transformation with 
large parameter p and majority rule, all of them defined on 2 X 2 blocks 9 . We find that 
the truncated Hamiltonian 7i' n is continuous at the transition line for the majority-rule 
transformation when restricted to a subspace containing a magnetic field and a nearest- 
neighbor interaction. On the other hand, we find that this is not the case for the deci- 
mation and large-p Kadanoff transformations restricted to the latter 2D subspace nor 
for the majority-rule transformation when restricted to the three-dimensional subspace 
containing magnetic field, nearest-neighbor and next-to-nearest-neighbor interactions. 
In all these latter cases, the renormalized magnetic field is non-zero, implying that the 
truncated RG map is discontinuous at the first-order transition line. Thus, the typi- 
cal situation seems to be that truncation induces discontinuities in the RG map when 
restricted to some finite-dimensional subspace of B 1 . 

This paper is organized as follows. In Section 2 we explain how the low-T expansions 
for renormalized observables can be obtained. We give three examples for the 2D Ising 
model: decimation, Kadanoff transformation and majority rule, all of them with block 
size 6 = 2. In Section 3 we study of these RGTs near the Ising first-order phase- 

8 In this case, the truncated Hamiltonian %' n = satisfies a fundamental property of the exact 

renormalized interaction W (if W exists at all), namely single-valuedness and continuity. 

9 It is known [11] that the first two transformations lead to non-Gibbsian renormalized measures at very 
low temperatures. This question is not clear for the last transformation, although it has been conjectured 
that it is also non-Gibbsian. Unfortunately, we are not aware of any real-space local transformation 
leading to a Gibbsian measure for the Ising model at low temperatures. In any case, we limit ourselves to 
studying question (A), for which the Gibbsianness or non-Gibbsianness is unlikely to play a significant 
role. 
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transition. Finally in Section 4 we present our conclusions. 



2 Series Expansions for Renormalized Operators 
2.1 Review of Low-T Expansions 

Let us consider for simplicity a ferromagnetic Ising model on a 2D square lattice. The 
spins take the values ±1 and interact through the Hamiltonian 

H = —K 5> i<7j . - 1) - H 5>,- - 1) (2.1) 

where the first sum is over all the nearest-neighbor pairs of spins, and the second one 
over every point i = (i X} i y ) of the lattice. The partition function for a system of N 
spins with periodic boundary conditions is then 



7 j N = ex p 

{a = ±l} 



K^ViVj-V + H^i-l) 



(2.2) 



We have absorbed the term f3 = 1/kT in the definition of the coupling constants K > 
and H. We are mainly interested in the zero-field case (H = 0), but for future 
convenience we keep the second term of the Hamiltonian (2.1). This term will be 
necessary to obtain the zero-field magnetization (see below). 

The first step to compute low-T expansions is to find out the ground states of the 
system at T = 0. In our case it is easy to realize that when H = there are only two 
translation-invariant ground states. Both of them are completely ordered configurations 
with magnetization +1 and —1 respectively. When H ^ then there is only one ground 
state whose magnetization is parallel to the magnetic field H . We will choose hereafter 
the (+l)-state as our ground state. This implies that the magnetic field should be 
always non-negative (H > 0). Furthermore, we have normalized the Hamiltonian (2.1) 
in such a way that 7Y(+1) = 0. 

Looking at eq. (2.2) it is easy to realize that each flipped spin is penalized by a factor 
A = exp(— 2H) in the partition function. And each unsatisfied bond (i.e. a bond with 
both spins in opposite states) is suppressed by a factor /i = exp(— 2K). All the spin 
configurations with n flipped spins and m unsatisfied bonds give the same contribution 
to the partition function (2.2) and equal to jj m X a . So we can group these configurations 
together and express the partition function as 

Z N ( f i,X) = J2Z { m N lfi m X n (2.3) 

where is the number of configurations with m unsatisfied bonds and n flipped spins 
that occur in the system. These numbers depend explicitly on the size of the system, 
as well as on the boundary conditions. The first term of the expansion corresponds 
to the ground state, the second to one flipped spin (n = 1, m = 4), the third to two 
nearest-neighbor flipped spins (n = 2, m = 6), and so on. With this choice of boundary 
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conditions, Z^\ = for odd values of m. This expansion is exact for finite N if all the 
2 N possible configurations are taken into account. 

The low-T expansion of the partition function (2.3) contains the most relevant terms 
when the temperature goes to zero. It can also be viewed as an enumeration of the low- 
energy excitations of the system. Here we are interested in developing an expansion valid 
as K — > oo with H bounded (i.e. an expansion in powers of /i (<C 1) whose coefficients 



m. 



are functions of A) 10 . Thus, the dominant terms are those with the smallest values of 
For a given value of m the possible values of n are finite. For excitations which do not 
see the boundary of the system the allowed values of n are given by n £ [m/4, m 2 /16] U 
[iV-m 2 /16,iV-m/4] (resp. [(m + 2)/4, (m 2 - 4)/16] U [N - (m 2 - 4)/16, N - (m + 2)/4]) 
when m/2 is even (resp. odd). All the terms with the same m, irrespective of n, are 
considered to contribute at the same order (i.e. A is considered to be of order 1). This 
feature implies that we can compute derivatives of the series expansions with respect to 
the magnetic field H . When the temperature is very close to zero only a few terms are 
needed to provide an accurate description of the system. 

Actually, the partition-function expansion is a technical tool to compute the ex- 
pectation values of some local operators: the energy density E = (c(o,o) cr (i,o)} an d the 
magnetization M = (c(o,o)}- The relations for a finite system are the following 

= 1 + ^it L = 1 -^iff = S^A" (2.4a, 

As before, the coefficients {Ef^^M^}} do depend on the lattice size and, in general, 
on the boundary conditions. 

Let us discuss now the thermodynamic limit (N — > oo) of these expansions. In 
this limit, the contribution of all the terms with the same m is not in general of the 
same order. In particular, for H > the configurations with n near N (for instance, 
n £ [N — m 2 /16, N — m/4] for m/2 even) are exponentially suppressed, and can therefore 
be dropped. Moreover, for H = the a — > —a symmetry implies that the contribution 
of the terms with n near zero is equal to the one of those with n near N. However, 
at H = + only the first set is selected. Therefore, for H > or H = + the correct 
expansion is obtained by taking all the terms with n near zero. 

On the other hand, the series corresponding to the partition function (2.3) are mean- 
ingless when N — > oo, as all the coefficients (except for Zq^ = 1) diverge in that 
limit. This is not true for the series (2. 4a, 2. 4b) whose coefficients have a well-defined 
limit 



£(/i,A) = ]T£ m , n /i m A n ; E m<n =]imEW (2.5a) 
M(/i,A) = ]TM m , n /i m A n ; M m , n = lim (2.5b) 



10 Different expansions are obtained when H —> oo and K remains bounded or when both K and H 
diverge with K/H — ► constant. 
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The limiting coefficients {E mjfl} M mjfl } do not depend on the boundary conditions of the 
finite systems. It is therefore reasonable to expect that {E mjfl} M mjfl } are the coefficients 
of the true infinite-volume series 11 . 

The series expansions for the zero-field case (H = + or A = 1~) can be easily 
obtained from the previous ones by summing over the index n. For example, M(/i) = 
£ m M ml j m where M m = £ n M m , a . 

In this paper we are mainly concerned about the computation of expectation values 
of more complicated local observables 0. By local operator we mean an operator which 
only depends on a finite number of spins. Our definitions of the energy density and the 
magnetization do satisfy this property. The previous procedure can be generalized to 
include also this case by adding to the Hamiltonian (2.1) a new term proportional to a 
translation-invariant version of the operator 0. 

However, this method is not feasible for very complicated operators, such as the ones 
considered in the next Section. In this paper we propose to use the definition 

(0) = Mml £ 0(a) e"^) (2.6) 
N ^°° Z N {a=±1} 

to overcome this problem. The term exp(— 7i) can be expanded in terms of configu- 
rations with m unsatisfied bonds and n flipped spins as we did in (2.3). In this case 
not all the configurations with the same values of m and n give the same contribution 
to the numerator of (2.6). This contribution is equal to jj m X a times the value of the 
operator 0(a) at the configuration. Let us consider a simple example. To compute the 
magnetization series one has to consider, for instance, the operator = C(o,o) (trans- 
lation invariance assures that the mean value of this operator will coincide with the 
magnetization (2.5b)). For instance, the contribution of the one-flip configurations is 
different depending on whether the flipped spin coincides or not with C(o,o)- I n the first 
case it is equal to —/i 4 X and in the second one to +/i 4 A. The same occurs for more 
complicated configurations (and operators). For a finite volume we obtain in this way 
an expansion similar to (2.4a, 2.4b). The final result (0) = £ m n O m ^ n jJ m \ n is obtained 
after performing the thermodynamic limit. 

The main advantage of this method is that it allows the computation of low-T series 
for arbitrary operators. Its main drawback is that we need to compute two series for 
each observable, not one as in the former method. 



2.2 Renormalization Group Transformations 

Let us begin by considering RGT as a map from an (infinite-volume) Gibbs measure v 
to a renormalized measure v' . Later on, the relationship between the measure v' and 
the renormalized Hamiltonian Ti 1 will be discussed. 

n A non-trivial interchange of limits is involved here, but it can presumably be justified rigorously 
by standard mathematical-physics techniques. We also need to be able to differentiate the free energy 
with respect to A at A = 1~, and this might be problematic since the free energy is not analytic in 
H at H = [14]. However, this function is infinitely differentiable, so we expect that our results are 
correct, as we compute everything in the stable phase. These issues arise in all uses of the conventional 
low-T-expansion technology, not only our own. 
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The first step is to define the renormalized spins. We divide the whole lattice into 
blocks (for simplicity we will assume here that these are 2x2 blocks). To each block 
Bi we associate a new (renormalized) spin a[. The RGT is the rule which gives the 
{a 1 } configuration from the original one {a}. This rule could be either deterministic or 
stochastic, but in any case the renormalized spin should only depend on the spins be- 
longing to the corresponding block (strict locality condition). Mathematically speaking 
we give a probability kernel T(cr, da 1 ). For each configuration of the original spins {c}, 
T(cr, •) is a probability distribution for the {a 1 } spins and furthermore, it satisfies the 
property J T(a } da 1 ) = 1. We assume that T is strictly local in position space and that 
it maps translation-invariant measures into translation-invariant ones. 

The probability distribution v' of the image system is therefore given by 

u' = uT = J dv(a)T(a 1 -) (2.7) 

and the expectation value of any local observable in this renormalized measure can be 
written as 

{0{a')) v , = I dv(a) \ ( T(a, da') 0(a')\ = {0(a)), (2.8) 



Thus, the probability kernel T when acting (to the left) on the measure dv(a) yields a 
probability distribution on the new spins {a 1 } (i.e. a renormalized measure v'). On the 
other hand, when T acts (to the right) on the operator O(a') the result is a composite 
operator 0(a) = (T ■ 0)(a) which depends only on the original spins. The expectation 
value of any local operator in the renormalized measure is equal to the mean value of a 
certain composite operator in the original measure. 

This discussion is general: the conclusions hold whether or not the systems can 
be described by a Hamiltonian 7i £ B 1 . Now we take into account the role of the 
Hamiltonians. Given an interaction 7i £ B 1 we can construct a measure on the spin 
configuration space using the Gibbs prescription 

dv(a) = dv°(a) \ e~ n ^ , (2.9) 

where dv°(a) is the a-priori measure we assign to the space of configurations of a 
single spin (in our case it is just the counting measure which gives to each state a 
probability 1/2). For finite systems, formula (2.9) gives the correct answer; but for 
infinite systems one has to be more careful and consider the limit of the finite-volume 
measures with given boundary conditions as the system size tends to infinity. For finite 
systems the relation between Hamiltonians and measures is one-to-one. However, in the 
thermodynamic limit this is not the case: one Hamiltonian can be associated to several 
measures (i.e. at first order phase transitions), and conversely there are perfectly sound 
measures which cannot be constructed via the Gibbs prescription from any Hamiltonian 

n £ b 1 [ii]. 

The Hamiltonian (2.1) does obviously belong to the set B 1 , so we can construct 
the measure v using (2.9). Then the expectation value (2.8) of any local renormalized 
operator can be written as 

(0)w = {6) v = lim J- £ 0(a) (2.10) 
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Here any given choice of boundary conditions gives rise to an original Gibbs measure v 
and a corresponding Gibbs measure v' . 

In Section 2.1 we showed how to obtain low-T expansions for a general mean value 
(0) v . Thus, the same procedure can be applied to (2.10), and series of the type (0) v > = 
J2m n 0'mnl jm ^ a are obtained. The practical applicability of this method relies heavily 
on the actual form of the kernel T as it is shown below. This procedure can also be 
easily generalized to several RG steps. 

It is important to remark that this method does not suffer from any of the pathologies 
which are exhibited by the RG when we try to define it as a map from a Hamiltonian 
space into a Hamiltonian space. Here we have not tried to define any renormalized 
interaction Ti 1 related with the renormalized measure v' via the Gibbs prescription (2.9). 
Our results are independent of the Gibbsian or non-Gibbsian nature of the renormalized 
measure. Let us illustrate this method with three examples: 



Example 1: Decimation 

This case is really simple because this transformation fixes one spin of the block to be 
the renormalized one. In particular, the (deterministic) kernel T takes the form 

T(a,a') = \{6(ala 2i ) (2.11) 

i 

where the product is over all sites i of the renormalized system. 

We are only interested in computing observables that are monomials of the spins 
(O = {c(o,o) 5 cr (o,o) cr (i,o)})- So it is enough to compute for each RGT the composite 
operator <7;. In this case this is equal to <7; = / T(<r, do')o[ = a 2 i- This implies that this 
case is trivial: the renormalized correlation functions are equal to the unrenormalized 
ones at twice the distance. These functions can be obtained in the 2D Ising model from 
the exact solution [15,16]. 



Example 2: Kadanoff Transformation 

This is given by the following (stochastic) probability kernel 



T(<t, a') = TT (2.12) 



where p is a free real parameter. Then, Oi = tanh [pJ2keB t a kj- The first terms can be 
computed by hand 

M'(/i,l") = tanh 4p - 4(tanh4p - tanh 2p)/i 4 - 4(3 tanh 4p - 2 tanh 2p)/i 6 

- (36tanh4p-4tanh2p)/i 8 + 0(/i 10 ) (2.13a) 
E'(n,l) = tanh 2 4p - 8(tanh 2 4p - tanh 4p tanh 2p)/i 4 

- 2(11 tanh 2 4p — 6 tanh 2ptanh 4p — tanh 2 2p)/i 6 

- (43 tanh 2 4p + 40 tanh 2p tanh 4p - 20 tanh 2 2p)/j s + 0(jj 10 ) (2.13b) 

The limit p — > corresponds to the case in which the a 1 are not correlated with the 
original spins and thus, the renormalized spins do not interact among them. For this 
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reason both quantities are zero. The limit p — > oo corresponds to the majority rule with 
equally-probable tie-breaker. This case will be treated in the next section. 

Example 3: Majority Rule 

In this case 

TV, a) = IJ S ^ - sign ^ £ a)j j (2.14) 

When sign(-) = we choose a' = — 1 or +1 with probabilities q £ [0,1] and 1 — q 
respectively. The composite operator a takes the form <t 4 - = sign (j2keB t a kj • The first 
terms for general q are: 

M'(/i,l-) = 1 - 8qfi 6 - (10 +Uq)fi 8 + 0{jj W ) (2.15a) 
E'(n,\) = 1 - 16g/i 6 - (20 + SSq - Aq 2 )/j 8 + 0(jj w ) (2.15b) 

The result with q = 1/2 was first reported in ref. [11]. Notice that the 0{jj a ) term 
vanishes. This is due to the fact that one-spin excitations cannot produce any flipped 
renormalized spin a' = —1. 



3 Study of the First— Order Phase Transition at Very 
Low Temperatures 

In MCRG calculations one chooses in advance a linear subspace V n C B 1 of the space 
of local Hamiltonians. Then, given certain renormalized expectation values, one tries 
to find a renormalized Hamiltonian 7i' n £ V n in such a way that a measure constructed 
from 7i' n is similar in some sense to the true renormalized measure v' . Most "recon- 
struction" methods are based on Schwinger-Dyson equations [17,18,19]. The idea is 
simple: minimize a certain functional (which depends on the method) involving both 
renormalized expectation values (the input) and renormalized couplings (the output). 
It can be shown [19] that these methods provide a unique solution 7i' n} which coincides 
with the true one Ti 1 if this latter interaction belongs to the trial subspace V n . The key 
property of these functionals is that they are strictly convex. 

Here we will consider the procedure given in ref. [11], Section 5.1.2. It is based on 
the minimization of the relative entropy density with respect to the true renormalized 
measure v' . This functional in also strictly convex and thus, the solution is unique in 
each V n if it exists. Van Enter et al. showed that the solution 7i' n should satisfy the 
following conditions 

{Oi) v , = {Oi) v , n VO t £ V n (3.1) 

where v' n is some Gibbs measure constructed from the Hamiltonian 7i' n . In this case we 
have the same number of equations as the number of unknown parameters. Note that 
7i' n might have multiple Gibbs measures; it is required that one of them satisfy (3.1). 

First-order phase transitions are characterized by the coexistence of several pure 
phases. Given one RGT, each pure phase is mapped to a renormalized measure 
v'^ k \ Given a subspace V n £ B 1 , we can solve eq. (3.1) for each renormalized measure 
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v'^ and obtain a corresponding truncated renormalized Hamiltonian Ti'^ . All these 
truncated Hamiltonians are uniquely defined and the important question is whether 
they are all equal or not. 

In the 2D Ising case, there are only two pure phases coexisting at T < T c 

and H = 0, and they are related by the a — > —a symmetry. For this reason, the 
even couplings in the truncated Hamiltonians are equal, and the odd couplings 

differ by a sign (in particular, the renormalized magnetic field in one phase is minus 
the renormalized magnetic field in the other one). Therefore, the truncated RG map 
is continuous and single-valued if and only if the renormalized Hamiltonians do not 
contain any odd interaction. Due to this symmetry, we only have to consider one phase 
(i.e. the (+l)-phase of Section 2 12 ). To check if the truncated Hamiltonian 7i' n has 
any odd term we simply solve (3.1) restricted to the even-coupling subspace of V n . If 
such a solution exists, then Ti'^ = Ti'^; if no such solution exists, then Ti'^ ^ 'H'^ ■ 



Case I. V 2 = {H, K} 

In this case 13 our subspace consists precisely of the Hamiltonians (2.1). We will try to 
match both the energy density and the zero-field magnetization by using a zero-field 
Hamiltonian at a different (lower) temperature. If this matching can be performed, it 
would mean that the truncated RGT for the subspace V2 is continuous at the transition 
line. 

Let us consider first the majority-rule map with random tie-breaker. We define K' 
as the nearest-neighbor coupling such that 

E'(K, 0) = E(K', 0) (3.2) 

Using the result (2.15b) and the well-known expansion of the Onsager solution 

£(/i,l) = 1 - 4/i 4 - 12/i 6 - 36/i 8 + 0{jj w ) (3.3a) 

M(/i,l") = 1-2//- 8/i 6 - 34/i 8 + 0(/i w ) (3.3b) 

we find that 

63 

= V2f + + 0(/) (3.4) 
The zero-field magnetization M at this particular temperature is equal to 

M(/i / ,l-) = l-4/i 6 -^/i 8 + 0(/i 9 ) (3.5) 

and this expansion should be compared with the renormalized magnetization M'(/i, A) 
given in (2.15a). We find that 

M(/j',r) > M'(/j,r) (3.6) 



12 To simplify the notation all the superscripts (+) will be dropped hereafter. 

13 One might consider [8] the even simpler case V\ = {H}. But this case is trivial since any RGT 
satisfying M(K, + ) 7^ M'(K, + ) is necessarily discontinuous at H = when restricted to the affine 
subspace A\ = V\ + (0, A') with fixed K ^> K c . 
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This equation means that we can account for the observed renormalized magnetization 
M'{n, 1") with a system at zero field and K' = -(1/2) log// 3K - (1/4) log 2. This 
system is not in a pure phase, but in a mixed phase because the renormalized mag- 
netization M'(/i,l~) lies strictly between ±M(//,1~). Thus, eq. (3.1) is satisfied by a 
measure v' 2 which is a convex combination of the two pure phases v^) characterizing the 
2D Ising model at inverse temperature K' and field H' = ± [i.e. v' 2 = av\^) + (1 — a)v^ K ) 
for a suitable a £ (0, 1)]. 

The same game can be played with the other two RGTs considered in Section 2. 
The easiest case is the decimation transformation, where conclusions can be drawn for 
every K > K c . In the 2D Ising model it is well-known that (c(o,o) cr (i,o)} > ( cr (o,o) cr (2,o)} 
for < K < oo. This implies immediately that E'(K,0) < E(K,0), and hence that 
K' < K if we take into account that E(K } 0) is a strictly increasing function of K . 
On the other hand, the renormalized magnetization coincides with the unrenormalized 
one (i.e. the RG flow follows the lines of constant magnetization). And M(K } + ) is 
also a strictly increasing function of K for K > K c . Combining both pieces we obtain 
that M(K',0+) < M'(K,0+) for all K > K c . This inequality is opposite to (3.6), 
because the direction of the RG flow is also opposite that of the majority-rule flow: 
it goes from low temperature to high temperature (K' < K). So, we have to increase 
the magnetic field to keep the magnetization constant, unless the magnetization at the 
starting point is zero. The latter condition only holds above the critical temperature. 
In summary, for any K > K c we cannot match the renormalized observables using a 
zero-field Hamiltonian. 

For the Kadanoff transformation and large (but finite) p the same result holds: one 
cannot match the energy densities and the magnetizations with a zero-field nearest- 
neighbor interaction. This can only be proved when p is large enough. The reason is 
clear: the leading term of E' is tanh 2 4p and if p is not large, then the solution of (3.2) 
does not satisfy //' <C 1 and the low-T series for // are then meaningless. 

For finite p we can always choose /i such that for /i < /i the leading term of E'(fi } 1) 
is dominated by a term which does not depend on /i. Then 

E'(fi, 1) = 1 - 4e- 8p + 0(e- 16p ) (3.7) 

if we choose /i ~ exp(— 3p). The solution of eq. (3.2) is then 

/i' = e~ 2p - h~ 6p + O(e- 10p ) (3.8) 

and 

M(fi', 1") = 1 - 2e~ 8p - 3e~ 12p + 0(e~ 16p ) (3.9) 

which should be compared with the expansion of the renormalized magnetization for p 
very large and /i < /i 

M'(fi,l-) = 1 - 2e~ 8p + 0(e- 16p ) (3.10) 

We find that the leading term of both quantities is the same, but the next-to-leading 
term is different. In particular we find that M'(/i, 1~) > M(/i', 1~), so we cannot match 
both E' and M' with a zero-field Ising interaction. This discussion is valid as long as p 
is large but finite. When p diverges the leading term of 1 — E'(fi } 1) is proportional to 
/i 6 and we re-obtain the result for the majority rule transformation with q = 1/2. 
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Case II. y 3 = {H, K, L} 

Now we are considering a Hamiltonian with an additional diagonal-next-to-nearest 
neighbor term L^ cr i cr fc- First of all we have to compute the renormalized mean value 
of the next-to-nearest neighbor correlation function. The result for the majority rule 
with random tie-breaker is 



F'(fji, 1) = (5(0,0)^(1,1)} = 1 - 4/i 6 - 64/i 8 - 336/i 10 - 1578/i 12 + 0(/i 14 ) 
We also need the functions H' and E' to the same order in /i. The result is 14 : 



(3.11) 



£(/i, 7 ,l) 
F(/i, 7 ,l) 
M(/i, 7 ,l-) 



M'(/i,l") = 1 - 4/i 6 - 32/i 8 - 168/i 10 - 816/i 12 + 0(/i 14 ) (3.12a) 
%1) = 1 - 8/i 6 - 63/i 8 - 312/i 10 - 1328/i 12 + 0(/i 14 ) (3.12b) 

The second step is to write down the expressions for (Oi) Vn , VO; £ V3. The result for 
zero magnetic field is 

1 _ 4/i Y - 12/i 6 7 8 - 24//A 12 - 32//A 10 + 36//A 8 
40/i 10 A 16 + 0(/i 8 A 6 ) (3.13a) 
1 _ 4/i Y - I6//7 8 - 36//A 12 - 40//A 10 + 36//A 8 
64/i 10 A 16 + 0(/i 8 A 6 ) (3.13b) 
1 _ 2/U 4 7 4 - 8/i 6 7 8 - 20//A 12 - 24//A 10 + 18//A 8 
- 40/i 10 A 16 + 0(/i 8 A 6 ) (3.13c) 

where 7 = exp(— 2L). Now we have to find out a pair (/i',7') such that 

£(//, 7 ', 1) = E\fi, 1); F(//, 7 ', 1) = 1) (3.14) 

The solution to leading order is // = 4/i 2 + £>(/i 4 ) and 7' = l/(32/i') 1/4 (l + £>(//)). This 
implies that K' « 2 A" - log 2 > and L' « (5/8) log 2 - K'/4 « (7/8) log 2 - K/2 < 0. 
So, as K — > 00, K' and — L' also diverge. The latter relation between /i' and 7' should 
be taken into account when computing the actual order of a given term in the expansion 
of the partition function Zjv(/u', 7', 1~) and its derivatives. In our case, this implies that 
the first two excitations to the ground state are of order /i' 3 and /i' 4 respectively. We 
have considered here all the excitations up to order 0(/i' 6 ). 

A straightforward computation leads to the next-to-leading terms 
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4/i 2 



1 



69 



-/i 2 + 72/i 3 + H|V + (9(/i 5 ) 



1 



1/4 



1 + 



327 



13/2 1 11 12 , /-/-)/ '5/2\ 



32/j' J f ' 256^ 16^2" ' 131072 
The magnetization (3.13c) computed at the latter solution is equal to 

2689 



M(/i', 7', 1") = 1 - 4/i 4 - 32/i 8 



16 



/i 10 + C(/i n ) < M'(/i,l") 



(3.15a) 
(3.15b) 

(3.16) 



This implies that we cannot match the renormalized expectation values with a zero-field 
interaction belonging to V3. 



14 We have computed this expansions by using a computer algorithm based on the Recursive Counting 
Method of ref . [20]. Further details can be obtained from the author. 
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4 Conclusions 



In this note we have shown how to compute low temperature expansions for the expec- 
tation values of local operators computed in the renormalized measure. In particular we 
have analyzed three RGTs: decimation, Kadanoff transformation with large but finite 
parameter p, and majority rule with random tie-breaker. All of them are defined on 
2x2 blocks. We have been able to compute the first terms of the series corresponding 
to the renormalized magnetization and nearest-neighbor two-point correlation function 
for all these transformations. 

The main goal of this note was the analysis of the truncation issue in the Ising model. 
The unrenormalized system is located at the Ising first-order transition line at very low 
temperature (H = 0, K ^> K c ). For the three transformations considered we have found 
that we need a magnetic field to solve the matching equations (3.1) when we restrict 
the truncated Hamiltonian to belong to a certain finite-dimensional subspace of B 1 . 
In particular, for the decimation and Kadanoff transformations this matching cannot 
be performed when restricting the equations to V2. For majority rule, in this case the 
equations admit a zero-field solution but when we consider the (larger) subspace V3 we 
also need a magnetic field. 

So its seems that truncation in the renormalized Hamiltonian induces some spurious 
odd operators (we have only found non-zero magnetic fields, but there is no reason 
why more complicated odd operators should not appear for larger subspaces V n ). So, 
these RGTs are discontinuous at the Ising transition line when restricted to some finite- 
dimensional subspace of the interaction space B 1 . 

However, these results do not clarify the interplay between truncation and non- 
Gibbsianness. It is known [11] that the decimation and Kadanoff transformations lead 
to non-Gibbsian renormalized measures when we start at low enough temperature; 
and in these cases we have shown that the truncated RGT are discontinuous. For the 
majority rule the situation is less clear, as the nature of the renormalized measure is not 
known. The authors of ref. [11] conjectured that in this case the renormalized measure 
is also non-Gibbsian, but they were able to prove it only for certain special block sizes 
(7 X 7, 41 X 41, . . .). In any case, this model leads to a continuous truncated RGT for 
the subspace V2, but a discontinuous one for V3. It is an open question what happens 
for larger subspaces V n . 

It would be very interesting to find a transformation which leads to a Gibbsian 
measure at low temperatures. In this case we could isolate the effect of truncation from 
non-Gibbsianness. A systematic study of the behavior of the estimates 7i' n is also an 
interesting problem, which deserves more attention in the future. 
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